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ABSTRACT 

Angular momentum transport in protostellar discs can take place either radially, 
through turbulence induced by the magnetorotational instability (MRI) , or vertically, 
through the torque exerted by a large-scale magnetic field that threads the disc. Us- 
ing semi-analytic and numerical results, we construct a model of steady-state discs 
that includes vertical transport by a centrifugally driven wind as well as MRI-induced 
turbulence. We present approximate criteria for the occurrence of either one of these 
mechanisms in an ambipolar diffusion-dominated disc. We derive "strong field" solu- 
tions in which the angular momentum transport is purely vertical and "weak field" 
solutions that are the stratified-disc analogues of the previously studied MRI channel 
modes; the latter are transformed into accretion solutions with predominantly ra- 
dial angular-momentum transport when we implement a turbulent-stress prescription 
based on published results of numerical simulations. We also analyze "intermediate 
field strength" solutions in which both modes of transport operate at the same radial 
location; we conclude, however, that significant spatial overlap of these two mecha- 
nisms is unlikely to occur in practice. To further advance this study, we have developed 
a general scheme that incorporates also the Hall and Ohm conductivity regimes in discs 
with a realistic ionization structure. 
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1 INTRODUCTION 

Tlie mechanisms responsible for transporting angular mo- 
mentum in protostellar discs are still not well understood. 
Perhaps the most generally relevant processes are: 

(a) Radial transport through turbulence induced by the 
magnetorotational instability (MRI; e.g. Balbus & Haw- 
ley 1998). The MRI acts, essentially, by converting the 
free energy of differential rotation into turbulent motions 
that transfer angular momentum radially outward via the 
Maxwell stress of small-scale, disordered magnetic fields. 

(b) Vertical transport by outflows driven centrifugally 
from the disc surfaces (e.g. Konigl & Pudritz 2000). The 
transfer of angular momentum from the inflowing gas to the 
wind is mediated by a large-scale, ordered magnetic field 
that threads the disc (likely corresponding to interstellar 
field lines advected inward by the accretion flow but possibly 
generated locally by a disc dynamo). 

Radial angular-momentum transport has been the hall- 
mark of generic viscous disc models, for which the MRI pro- 
cess is thought to provide the most cogent physical basis. 
On the other hand, the ubiquity of bipolar outflows and jets 
in protostars and the apparent correlation between outflow 
and accretion signatures in these systems (e.g. Hartigan et 
al. 1995) have led to the suggestion that disc outflows might 



be a key ingredient of the accretion process on account of the 
high angular-momentum transport efficiency of centrifugally 
driven winds (e.g. Konigl 1989). The possibility of a central 
role for vertical angular-momentum transport is supported 
by recent HST observations of protostellar jets, which have 
been interpreted in terms of magnetocentrifugal winds car- 
rying a large fraction of the excess angular momentum of 
the accreted matter from the inner regions (r < 1 AU) of 
the associated discs (see Ray et al. 2006 for a review). Al- 
though thermal pressure may contribute to the driving of the 
outflows (e.g. Pesenti et al. 2004), the stresses exerted by a 
large-scale magnetic field dominate the flow acceleration and 
coUimation (Blandford & Payne 1982, hereafter BP82; see, 
however, Soker & Regev 2003 for an alternative view of the 
importance of thermal acceleration) and are the key to the 
efficient vertical angular-momentum transport indicated by 
the HST measurements. 

Previous analytic and numerical investigations have 
considered the above two mechanisms under various sim- 
plifying assumptions. MRI studies have examined both the 
linear and nonlinear stages of the instability (e.g. Sano & 
Stone 2002a,b; Salmeron & Wardle 2003, 2005). The long- 
term global evolution of the instability in vertically strat- 
ified discs that are threaded by a single-polarity magnetic 
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field is still an open question (e.g. Miller & Stone 2000). 
The study of wind-driving discs has been complicated by 
the need to solve simultaneously for the inflow and outflow 
structures and by the non-local nature of the MHD solution. 
Progress has nevertheless been made semi-analytically by 
adopting a self-similarity formulation (e.g. Wardle & Konigl 
1993, hereafter WK93; Li 1995, 1996; Ferreira 1997) as weU 
as numerically, with recent simulations reporting an evolu- 
tion to a quasi-stationary state (e.g. Casse & Keppens 2002; 
Kuwabara et al. 2005). Since protostellar discs are weakly 
ionized over most of their extent, magnetic diffusivity ef- 
fects have to be taken into account in modeling the accre- 
tion flow in both cases. In particular, both the MRI-driven 
turbulence and the centrifugal wind-launching mechanism 
require a minimum level of field-matter coupling to be ef- 
fective (e.g. WK93; Sano & Inutsuka 2001). Therefore, the 
detailed ionization structure and conductivity properties of 
the disc need to be determined in order to incorporate either 
one of these processes into a realistic model. 

Angular momentum removal from real discs is likely to 
involve both the radial and vertical transport mechanisms 
discussed aboveQ There have been a few attempts in the lit- 
erature to construct quasi-steady disc models in which both 
of these mechanisms operate (e.g. Lovelace et al. 1994; Casse 
& Ferreira 2000; Ogilvie & Livio 2001), but these treatments 
invariably used a prescription for radial angular-momentum 
transport that did not explicitly account for its origin in 
MRLinduced turbulence. However, the fact that disc regions 
in which the entire excess angular momentum is carried away 
by a centrifugally driven wind are generally stable to the 
fastest growing mode of the MRI (WK93) indicates that 
a more refined approach is required. Broadly speaking, the 
MRI is suppressed in regions where the ratio a = vaz/cs of 
the Alfven speed corresponding to the vertical field compo- 
nent Bz and the isothermal sound speed is comparatively 
large (> 1); it may be expected to evolve into a fully devel- 
oped turbulence when a <^ 1 (provided also that the gas is 
well coupled to the field) . Since the value of a increases with 
distance z from the disc midplane, these considerations sug- 
gest that it could in principle be possible for both vertical 
and radial angular-momentum transport to take place at the 
same disc radius r, with the radial transport confined to a 
given range in z and with the vertical transport dominating 
at greater heights. 

The purpose of this paper is to construct semi-analytic 
models of protostellar discs that incorporate both centrifu- 
gally driven winds and MRI-induced turbulence and to iden- 
tify the relevant parameter regimes where either one of these 
two mechanisms dominates the angular momentum trans- 
port. We pay particular attention to the possibility that 
both of these mechanisms operate at the same disc radius 
and consider how they might affect each other in that case. 
We describe our wind-driving disc model (formulated fol- 
lowing WK93) and present an illustrative solution of a disc 
with a predominantly vertical angular momentum transport 



^ Vertical angular momentum transport could alternatively take 
place through the process of magnetic braking, which involves 
the propagation of torsional Alfven waves into the surrounding 
interstellar medium. This possibility is neglected in what follows, 
but see Krasnopolsky & Konigl (2002). 



in Section 2. In Section 3 we outline our method of including 
MRI-induced turbulence (taking account of recent numeri- 
cal simulations of the nonlinear evolution of the MRI) and 
illustrate it with a solution of a disc with a predominantly 
radial angular-momentum transport. We then consider the 
case of combined transport by wind and turbulence in Sec- 
tion 4 and summarize our results in Section 5. 



2 VERTICAL TRANSPORT 

We have devised a general scheme for modeling wind-driving 
protostellar discs using a tensor-conductivity formalism and 
a realistic vertical ionization structure0 The conductivity 
tensor makes it possible to include in a systematic way the 
three relevant magnetic diffusion mechanisms (ambipolar. 
Hall and Ohmic). The solutions presented in this paper are, 
however, derived under the simplifying assumptions of pure 
ambipolar diffusivity and (except in Section 4) constant ion 
density. This enables us to apply the analytic results ob- 
tained for this case in WK93 to the formulation of a pre- 
scription for MRI-induced turbulent transport in a wind- 
driving disc model. A wider class of solutions of our general 
scheme will be presented in future publications. 

Our wind solutions are based on the magnetocentrifu- 
gal mechanism for "cold" disc outfiows proposed by BP82. 
In this scenario, matter at the disc surface is accelerated by 
magnetic stresses if the open magnetic field lines that thread 
the disc are inclined at a sufficiently large angle to the rota- 
tion axis (the "bead-on-a-wire" effect), potentially reaching 
super- Alfvenic speeds. Our modeling procedure is based on 
that of WK93 and we refer the reader to that study for 
further details. In brief, we solve the mass and momentum 
conservation equations for the neutral gas, and the induction 
equation for the evolution of the magnetic field, assuming a 
steady-state, geometrically thin and nearly Keplerian accre- 
tion disc and neglecting all radial derivatives except that of 
W0. The disc is further assumed to be isothermal and weakly 
ionized. The latter condition implies that ionization and re- 
combination processes, as well as the inertia and thermal 
pressure of the charged particles, have a negligible dynam- 
ical effect. In this limit, separate equations for the ionized 
species are not needed, and their effect on the neutrals can be 
conveyed by means of a conductivity tensor cr(r, z) (assum- 
ing cylindrical polar coordinates; see Wardle 1999 and refer- 
ences therein). In addition, the field must satisfy V ■ B — 
(which in the thin-disc approximation implies that Bz can 
be assumed not to vary with height in the disc), whereas the 
current density J must satisfy Ampere's and Ohm's laws. 

The foregoing equations are integrated vertically from 
the midplane (indicated by a subscript 0) up to the position 
Zs of the critical (sonic) point, and the values of the physical 
variables there are estimated. The solution is then integrated 

^ The main ionization sources outside the innermost ~ 0.1 AU 
from the central object are nonthermal: interstellar cosmic rays, 
X-ray and UV radiation emitted by the magnetically active pro- 
tostar and the decay of radioactive elements present in the disc. 
The gas near the disc surface is ionized primarily by X-rays in 
the vicinity of the star, and by UV photons (with a contribution 
also from the interstellar radiation field) in the outer disc (e.g. 
Semenov et al. 2004; Glassgold et al. 2005). 
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backward to a fitting point and iterated until full conver- 
gence is reached. The resulting disc solution is matched on 
to a global BP82-type self-similar wind solution, obtained 
by imposing the regularity condition at the Alfven critical 
point. This final step effectively determines the value of 
at the disc surface (subscript b, defined as the height where 
equals the Keplerian speed vk)- The overall solution is 
determined by the following parameters: 

(i) ao = vaz,o/cs, the midplane ratio of the Alfven speed 
(based on the uniform vertical field component) to the sound 
speed, which measures the magnetic field strength. 

(ii) 77, the ratio of the Keplerian rotation time f2^^ = 
r/vK to the neutral-ion momentum-exchange time (see 
WK93); this parameter measures the degree of coupling be- 
tween the neutrals and the ionized gas component (and mag- 
netic field) and does not vary with height if the ion density 
is const ant 

(iii) e = —Vro/cs, the normalized inward radial speed at 
the midplane. In practice, the value of this parameter is 
determined in our solutions (for given values of the other 
parameters) by imposing the Alfven critical-point constraint 
on the wind. 

(iv) Cs/vk = h^/r, the ratio of the tidal scale height to 
the radius - a measure of the geometric thinness of the disc. 

(v) €b = —cE^/csBz, the normalized azimuthal compo- 
nent of the electric field E. This parameter vanishes in a 
strictly steady-state solution but is nonzero if the magnetic 
field lines drift radially. Our formulation remains valid if the 
field-line drift speed is small compared to \vr\ (see WK93). 

Q 

Fig. [1] shows an example of a local disc solution that 
matches on to a global BP82 wind solution. This solution 
resembles the one presented in Fig. 10 of WK93, although 
it was derived for a different choice of the parameters 77 and 
ao. Both of these solutions illustrate the main properties 
of wind-driving, strongly magnetized (ao < 1) discs. In the 
quasi-hydrostatic layer just above the midplane the mag- 
netic field lines are (radially) bent and (azimuthally) sheared 
and the field removes angular momentum from the matter. 
The scale height h in this region is reduced by magnetic 
squeezing, which dominates the tidal compression. Above 
this layer lies a transition region where the magnetic en- 
ergy comes to exceed the gas internal energy and the field 
lines are locally straight. In both of these zones, the fluid 

^ In the more general formulation, one needs to specify the values 
of the independent components of the conductivity tensor - viz. 
the field-aligned (en), Pedorsen (crp) and Hall (cth) conductivi- 
ties. In this case, the magnetic coupling is expected to decrease 
with height in the ambipolar diffusion-dominated upper section 
of the disc, where the ion density is most strongly affected by 
the decrease of the neutral gas density with z (e.g. Salmeron & 
Wardle 2005). In the simplified approach adopted in this paper 
we use the parameter rj = (47rCg /c^)a^(crp/f2K) instead (with 
a oc p"^/^ increasing with 2) and (except in Section 4) assume 
that it remains constant with height. 

^ Our ansatz 63 = effectively fixes the value of Br at the disc 
surface. In a generalization to a global treatment, B^^h would be 
determined by the radial distribution of B^ along the disc (e.g. 
Ogilvie & Livio 2001; Krasnopolsky & Konigl 2002), and the value 
of eB could be inferred from the overall magnetic field structure 
and the disc conductivity. 




0.5 1 1.5 2 2.5 
z/h^ 

Figure 1. Vertical structure of a strongly magnetized, wind- 
driving disc. Top: Br, \B^\, and p. Bottom: velocity components. 
The curves terminate at the sonic point. The disc solution param- 
eters are indicated in the figure. The corresponding parameters 
of the BP82 self-similar wind solution are k = 3.2 X 10~* (nor- 
malized mass-to-magnetic flux ratio), A = 395 (normalized total 
specific angular momentum), and = = 1.46. 



azimuthal speed is sub-Keplerian. Since the gas angular ve- 
locity decreases with radius, the field lines (which move with 
a constant angular velocity) eventually overtake the mat- 
ter and fling it out centrifugally. The region where Vr and 
{v^ — Vk) become positive constitutes the base of the wind. 
It has been shown (WK93; Konigl 1997) that solutions of 
this type satisfy 

(277)-'/' <ao<V3<er,< vk/2c, . (1) 

The first inequality in equation ([ij expresses the require- 
ment that the fluid remain sub-Keplerian below the wind- 
launching region, the second gives the wind launching con- 
dition {Br,b/Bz > 1/^3), and the third gives the constraint 
that the base of the wind be located above a density scale 
height. The first two conditions imply a minimum value for 
the field-matter coupling parameter i] (which a more de- 
tailed analysis sharpens to > 1), whereas the last two 
imply that magnetic squeezing dominates the vertical con- 
finement of the disc {h < /it). The last inequality specifies 
that the rate of ambipolar diffusion heating at the midplane 
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Figure 2. Structure of a weakly magnetized disc with verti- 
cal magnetic angular-momentum transport in the limit where no 
wind develops (channel solution) . The dashed line marks the nom- 
inal boundary of the MRI- unstable zone discussed in the text. 



should not exceed the rate of release of gravitational poten- 
tial energy. 



3 RADIAL TRANSPORT 

Numerical simulations indicate that vigorous radial angular 
momentum transport is sustained by the MRI when the ini- 
tial (subscript i) Elsasser number Ai = i/»?ohm^^K is -SJ 1 
(e.g. Sano & Inutsuka 2001; Sano & Stone 2002a,b), where 
VAz,i is the (midplane) Alfven speed based on the initial ver- 
tical field Bz,i and 77ohm is the Ohmic diffusivitylfl Now, for 
the vertical initial field geometry assumed here, the currents 
and fields are mutually orthogonal and ambipolar diffusion 
acts as a field-dependent Ohmic resistivity (e.g. Balbus & 
Terquem 2001). Under these conditions one can show that 
Ai = rji, so the minimum-coupling condition in this case 
is the same for vertical and for radial angular-momentum 
transport. 

What determines which mode of angular momentum 
transport dominates when the coupling is sufficiently strong 
{v si 1)? As already noted in Section 1, the MRI is sup- 
pressed for sufficiently high values of the field-strength pa- 
rameter a. A useful clue to the lower limit for MRI stability 
is provided by the first inequality in equation ([l]). As dis- 
cussed in WK93, when this inequality is violated the surface 
layers of the disc become super-Keplerian, implying outward 
streaming motion that is unphysical in the context of a pure 
wind-driving disc model. However, as elaborated on below, 
an outward motion of part of the disc material is a charac- 
teristic of the two-channel MRI mode that underlies MRI- 
induced turbulence. Further support for this choice of an 
upper limit on a in the MRI-unstable region is provided by 
the fact that wind-driving disc solutions that satisfy the con- 
straints expressed by equation ^ are stable to the fastest 
growing linear mode of the MRI, essentially because the min- 
imum wavelength of this mode exceeds the magnetically re- 



^ A is often referred to in the MRI literature as Rbm, the mag- 
netic Reynolds number, but in general Rcm ^ A. 



duced disc scale height (see WK93)0 The above limit on 
a was applied by WK93 at the midplane, since they con- 
sidered discs in which the entire excess angular momentum 
of the accretion flow is transported by a wind. As we are 
now concerned with the possible existence of an alternative 
mode of transport, we propose to generalize this condition 
and apply it locally (i.e. using a{z) oc p{z)~^^^ instead of its 
midplane value) to differentiate the MRI-unstable section of 
the disc (where 2ria^ < 1) from the region where only verti- 
cal angular momentum transport takes place (2r;a^ > 1). 

By integrating the disc equations of the model described 
in Section 2 for the case in which a < {2rj)~^^^ over the bulk 
of the disc we have verified that such systems do not develop 
winds. Since our current formulation does not allow for an al- 
ternative means of transporting angular momentum through 
the disc surfaces (such as magnetic braking, "failed" winds 
that do not cross the relevant critical surfaces, or nonsteady 
phenomena), the only remaining option is for the angular 
momentum to be transported magnetically between differ- 
ent heights and then to be carried away radially by bulk 
fiow (with the gas that loses angular momentum moving 
in and the material that gains angular momentum moving 
out). Fig. [2] shows a solution of this type. The radial field 
component oscillates in the region where 2ria^ < 1, with 
exhibiting corresponding oscillations with an opposite sign. 
The regions where BzdB^/dz < lose angular momentum 
to the field and have Vr < 0, with the converse behaviour 
characterizing regions where BzdB^/dz > 0. We associate 
this solution with a two-channel MRI mode for the case of 
a stratified, non- ideal MHD discH] As was first shown by 
Goodman & Xu (1994), such modes (which they considered 
in the unstratified, ideal-MHD case) are exact solutions in 
both the linear and the nonlinear regimes. Their analysis, 
however, revealed that these solutions are unstable to para- 
sitic modes, and subsequent numerical simulations (carried 
out in both the ideal and non-ideal regimes) have verified 
that they rapidly evolve into a turbulent state (e.g. Hawley 
et al. 1995; Sano et al. 2004, hereafter S04). 

To account for the turbulent angular momentum trans- 
port that we expect to develop in disc regions where ri > 1 
and 2a^rj < 1, we adopt the results of numerical simula- 
tions, which have found that the main contribution to the 
r(j> component of the stress tensor is provided by the tur- 
bulent Maxwell stress and that its spatial and temporal av- 
erage <^ Wrcf, ^ can be expressed in terms of the similarly 
averaged magnetic pressure by 

<w.0>^O.5 ^f (2) 
on 

(e.g. Hawley et al. 1995; S04). Using the characteristic ratios 

< B| > / < > « 24 and < ^ / < > « 3 for 
MRI turbulence reported in Table 4 of S04, we substitute 

< > 28 < > into equation © to get 

B^ 

<Wr^>^i4y-^^ , (3) 

where F = < Bf > /B^i = Y{f3) is listed in Table 2 of S04 

^ Note, however, that the formal demonstration of this fact re- 
quires using only the third inequality in equation JTJ. 
^ Similar field configurations were obtained by Ogilvie & Livio 
(2001), who also related them to the channel solution. 
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Figure 3. Structure of a weakly magnetized disc with purely 
radial angular-momentum transport given by the MRI-induccd 
turbulence prescription discussed in the text. The curves termi- 
nate at the point where the parameter a satisfies 2ria? = 1. 

for a range of isothermal, uniform-Bz,i models as a function 
of /3i = , the initial plasma (3 parameter. As the values 
of 13 in our applications are smaller (by up to ~ 1 — 2 orders 
of magnitude at the midplane) than the ones listed in S04, 
we resort to extrapolation, which can be justified by the 
relatively weak dependence of Y on /3: we obtain Y = 2, 1 
and 10 for jS — 10, 100 and 500, respectively. 

Radial angular-momentum transport is incorporated 
into the disc angular momentum conservation equation by 
adding (in the relevant disc region) the turbulent stress 
term, which we do by approximating (l/r)9(r^ ^ Wr^, 2> 
)/dr Wr,/, 3>. This equation then reads 

-^ + P''^^ = —c r ' 

where we used d{rv^)/dr ~ vk/2 and « 0. We as- 
sume that the small-scale (turbulent) magnetic field com- 
ponents only contribute to the radial angular-momentum 
transport. We interpret all the other magnetic-field (and ve- 
locity) terms that appear in the disc structure equations as 
corresponding to the mean values of the respective quan- 
tities. In a similar vein, we interpret the first term on the 
r.h.s. of equation Q as representing the vertical angular- 
momentum transport induced by the mean field. This rep- 
resentation is akin to that of previous models that in- 
corporated radial angular-momentum transport into discs 
threaded by an open magnetic field (see Section 1). The 
underlying assumptions of this approach - that the rela- 
tionships obtained from local shear-box simulations remain 
apphcable in a stratified disc, that <^ BrBz 3> averages out 
to zero and that ^ {B^ + B^) ^ does not have a significant 
z dependence - will, however, need to be verified by means 
of global numerical simulations. According to the above for- 
mulation, vertical angular-momentum transport associated 
with a large-scale, ordered field occurs throughout the disc 
(with the angular momentum deposited into the mean field 
eventually transferred back to matter in the wind region 
above the disc), whereas radial transport effected by MRI- 
induced turbulence is confined to the layer where the condi- 
tions 1 < ?? < (2a^)"^ are satisfied. 

Fig. Oshows a weak-field disc solution, obtained by im- 
plementing the above prescription for the same values of -q 



and ao as in the channel solution presented in Fig. (2] In 
this case the angular momentum transport is purely radial 
throughout the disc. The azimuthal field component at the 
top of the disc does not generally vanish in weak field disc 
configurations of these type, implying a net vertical torque 
(oc BzB,^,b) on the disc. However, in the current formula- 
tion such a torque can only be included if the wind launch- 
ing condition is satisfied. Since this is not the case here, we 
chose to eliminate all the magnetic terms from the derived 
solution for self consistency. The resulting error is, however, 
small since the mean-field contribution to the disc dynamics 
(in particular, to angular momentum transport and to ver- 
tical squeezing) is rather minor in this instance. The lack of 
strong magnetic squeezing is refiected in the more gradual 
decline of the density with height in this solution as com- 
pared with the pure wind-driving disc model presented in 
Fig.IIl 



4 COMBINED TRANSPORT 

We now consider the case of a disc with moderate coupling 
(77 ~ 1) and field strength (ao < (277)"^'^^ but a ~ 1 at 
the surface), so both radial angular-momentum transport 
by a turbulent magnetic field (between the midplane and 
the height zt where 2rja? ~ 1) and vertical transport by 
a large-scale field (between z = Q and the top of the disc 
at Zh) take place at the same radial location. Solutions of 
this type occupy a limited region of parameter space, cor- 
responding to a large enough ao for a wind to be launched 
from the disc surfaces and a sufficiently small rjo to render 
2r7oao < 1 (but 770 > 1 so as not to suppress the growth of 
MRI-unstable modes). We have found that if r; is assumed 
to be both > 1 and constant with height (as was done in the 
previous Sections) , the vertical extent Azmri over which ra- 
dial transport operates is very narrow. The value of Azmri is 
limited by the fact that ao has a lower limit below which the 
wind-driving surface layers of the disc become sub-Keplerian 
(corresponding to an effective violation of the first inequality 
in equation [T] in these regions) . To circumvent this limita- 
tion, we have allowed to fall slightly below 1 (with the 
expectation that it is still large enough to sustain MRI tur- 
bulence; see Sano & Stone 2002a,b) and to decrease with 
height. The adopted profile is motivated by the inferred be- 
haviour of 77(2) in the ambipolar diffusion-dominated upper 
layers of realistic protostellar discs at r ~ 1 — 10 AU, where 
the ion density is most strongly affected by the decrease of 
the neutral gas density with height (see footnote |3]) . 

Fig. |4] shows illustrative solutions obtained under these 
assumptions. The solid lines depict the structure of the disc 
when only vertical angular-momentum transport is included 
(cf Fig. [T|. The dashed lines indicate how the above so- 
lution is modified when radial transport is incorporated 
using the prescription discussed in Section 3 and keeping 
the value of ao fixed. Key properties of the two solutions 
are compared in Table [TJ assuming that Cs and po also re- 
main unchanged. The mass accretion rate per disc circum- 
ference is evaluated from Mi-a/2-nro — —2 J^^ pVrdz; the 

listed value (A4in) is this quantity normalized by (poCa/iT). 
The vertical and radial torques (per unit area) on the disc 
are normalized by toBq/A-k and are given by Tz = kA 
and Tr = (47r<c: Wr^, 2>/Bo)(A2mri/''o), respectively (with 
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Table 1 . Comparison of the two solutions shown in Fig. |4] 



Figure 4. Wind-driving disc with a moderately strong field (ao = 
0.55) and rj decreasing from 0.63 to 0.5 between the midplane 
and the surface. Density, magnetic field components (top) and 
velocity components (bottom) are shown for a solution that only 
incorporates vertical angular-momentum transport (solid lines) 
and also for a case where radial transport is included (dashed 
lines) in the region where 2ria^ < 1 (to the left of the vertical 
dotted line). The curves terminate at the respective sonic points. 
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k) and a correspondingly lower vertical torque (measured by 

The results shown in Fig.|4]illustrate the possibility that 
the incorporation of MRI-induced turbulence could increase 
the overall angular momentum transport in the disc, but re- 
duce the amount of angular momentum carried away by the 
wind. It may be worth noting in this connection that the 
enhanced squeezing of the disc brought about by the addi- 
tion of radial transport also reduces the width of the tur- 
bulent layer in comparison with the value of zt that would 
have been estimated from the pure wind solution (in which 
a increases less rapidly with z). These considerations sug- 
gest that the joint operation of the two angular-momentum 
transport mechanisms might be self-limiting, in which case 
it would be unlikely to introduce a dynamical instability. 
This conclusion could be tested through an explicit stabil- 
ity analysis (e.g. Konigl 2004 and references therein) or by 
means of numerical simulations. 



A^MRi = Zt for the given setup), whereas the azimuthal field 
component at the disc surface is obtained from |-B,^,bj/Bz = 
k(A — 1). The heights zt, Zh, and Za are listed in units of the 
scale height h, which is determined for each solution as the 
value of z where the density drops to po/\/e- 

The addition of radial angular-momentum transport re- 
sults in this case in an increased inflow speed (measured by 
e) and a higher mass inflow rate (Min)- The higher \vr\ cor- 
responds to a stronger radial neutral-ion drag (since the 
ions and fleld remain pinned at the midplane under the as- 
sumption that eB ~ 0; see footnote |4)| , and therefore (by the 
ion radial force-balance condition; see WK93) to a stronger 
bending of the magnetic field lines away from the vertical 
(reflected in the higher value of Br,h)- This, in turn, increases 
the level of magnetic squeezing (reflected in the lower value 
of h), resulting in a stronger density stratification. Although 
the thickness of the disc remains virtually unchanged (at 
Zb = 1.25 /it) and the height of the sonic point actually goes 
down in the combined solution (from 2.90 /iT to 2.82 /it), 
the increased stratification (measured by Zh/h and Zs/h, re- 
spectively) means that the density at the top of the disc is a 
lower fraction of the midplane density than in the pure wind 
solution, leading to a lower wind outfiow rate (measured by 



5 CONCLUSION 

We have constructed a semi-analytic scheme for modeling 
magnetized protostellar accretion discs. The model takes 
into account both vertical angular-momentum transport as- 
sociated with a large-scale magnetic fleld that threads the 
disc and drives a centrifugal wind from its surfaces and radial 
angular-momentum transport through MRI-induced turbu- 
lence. A novel feature of this model is an explicit prescrip- 
tion that identifies the vertical extent of the disc region that 
is susceptible to sustainable MRI turbulence at any given 
radius and incorporates a turbulent-stress term (based on 
published results of numerical simulations and related to the 
local magnitude of the large-scale fleld) into the angular- 
momentum conservation relation for this region. Our gen- 
eral scheme employs a tensor-conductivity formalism and 
incorporates a realistic vertical ionization structure for the 
disc. However, in this paper we specialize to the ambipolar 
diffusion-dominated regime and assume that the ion density 
is nearly constant with height. This allows us to compare our 
results with the pure vertical-transport solutions obtained 
for this case by WK93 and to use the algebraic relations de- 
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rived in that paper in the analysis of the combined vertical 
and radial transport. 

We have concentrated on the case where the neutral- 
ion coupling parameter ?? is > 1, which allows both the 
wind-driving and MRI-turbulence mechanisms to operate 
efficiently. Vertical transport through a wind is expected 
to dominate when the field-strength parameter a (the ra- 
tio of the local Alfven and sound speeds) is sufficiently high 
(< f), whereas radial MRI-induced transport should prevail 
for a <^ 1. We constructed a low-a disc solution exhibit- 
ing oscillatory field and velocity profiles and suggested that 
it evidently represents the steady-state analogue in a diffu- 
sive (and stratified) disc of the nonlinear MRI channel mode 
originally identified by Goodman & Xu (1994). Although 
the question of the stability of stratified channel solutions 
in the non-ideal MHD regime has not yet been explicitly 
addressed, the current evidence (from linear analyses and 
numerical simulations) suggests that they may well be un- 
stable. On the assumption that, in fact, such flows rapidly 
become turbulent, we incorporated a turbulent-stress pre- 
scription (based on the results of numerical simulations) into 
the angular momentum conservation equation for this re- 
gion. We showed that, when implemented, this prescription 
results in an accretion-disc solution featuring pure inflow. 
We argued that radial angular-momentum transport termi- 
nates when 2r)a? increases above ~ 1, and we presented a 
hybrid solution in which this transition occurs at a height zt 
below the disc surface (^b), with radial transport (restricted 
to the layer z < Zt) augmenting the vertical transport (which 
operates for all z < Zb)- We found that the joint operation of 
these two mechanisms at the same radial location in the disc 
is apparently self-limiting. In particular, we showed that the 
incorporation of radial transport, while increasing the total 
amount of angular momentum removed from the inflowing 
gas, leads to a reduction in the angular momentum flux car- 
ried away by the wind. 

Our results strongly indicate that the two angular- 
momentum transport mechanisms are unlikely to have a sig- 
nificant spatial overlap under realistic circumstances. This 
is because of the stringent constraints that must be satis- 
fied for this to occur: the parameter a must be small enough 
(< {2rj)~^^^) near the midplane to allow MRI turbulence 
to develop over a significant vertical extent, but it cannot 
be so small that the field in too weak to launch a wind (or 
else that the surface layers become sub-Keplerian) . In prac- 
tice, these constraints require moderate midplane values of 
{o-o ^ 0.5) and rather low (and slowly decreasing with 
height) values of ri {rio < 1), which are unlikely to be at- 
tained over measurable radial extents in real discs. If this 
result indeed applies generally to protostellar discs in the 
ambipolar diffusion regime then it is of interest to inquire 
where each of these mechanisms is likely to dominate in such 
systemslfl As an example, we use the normalizations of ao 
and -qo given (on a scale of ~ 100 AU) by equations (3.50) 
and (3.51), respectively, of WK93 and the radial scalings 
{B oc r"^''*, p oc r~^, ion density pion oc p^^^) correspond- 



° It is, however, worth pointing out again in this connection 
that vertical angular-momentum transport mediated by large- 
scale fields could take place (e.g. in the form of magnetic braking) 
even if the conditions for wind launching are not fulfilled. 



ing to the asymptotic rotational core-collapse solution of 
Krasnopolsky & Konigl (2002), which imply that 2?7oao cx r 
drops below 1 at r ~ 40 AU (with disc winds dominating 
the angular momentum transport at larger radii and MRI- 
induced turbulence turning on closer in). 

The above example should, of course, be taken only as 
an illustration, since the radial parameter profiles in actual 
discs may be quite different. Furthermore, as the density 
and column density in the disc increase with decreasing ra- 
dius, the Hall conductivity regime, and eventually also the 
Ohmic regime, will be encountered (e.g. Balbus & Terquem 
2001), with different regimes expected to dominate at differ- 
ent heights at any given radial location (Wardle 2004). We 
are currently in the process of applying our general scheme 
to the investigation of the various conductivity regimes in re- 
alistic protostellar disc models. The results will be reported 
in future publications. 
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